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ABSTRACT 

Motivated by recent detection of transiting high-density super-Earths, we explore the detectability 
of hot rocky super-Earths orbiting very close to their host stars. In the environment hot enough for 
their rocky surfaces to be molten, they would have the atmosphere composed of gas species from the 
magma oceans. In this study, we investigate the radiative properties of the atmosphere that is in the 
gas/melt equilibrium with the underlying magma ocean. Our equilibrium calculations yield Na, K, 

Fe, Si, SiO, O, and O 2 as the major atmospheric species. We compile the radiative-absorption line 
data of those species available in literature, and calculate their absorption opacities in the wavelength 
region of 0.1-100 /im. Using them, we integrate the thermal structure of the atmosphere. Then, 
we find that thermal inversion occurs in the atmosphere because of the UV absorption by SiO. In 
addition, we calculate the ratio of the planetary to stellar emission fluxes during secondary eclipse, 
and find prominent emission features induced by SiO at 4 /im detectable by Spitzer, and those at 10 
and 100 /mi detectable by near-future space telescopes. 
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1. INTRODUCTION 

We have recently entered a new era of characterizing 
the atmospheres of low-mass exoplanets with masses of 
< 30 Mg. Transmission and/or emission spectra have 
been ob tained for several low-mass exoplanet s in tran¬ 
sit fe.g.. iDeming et all 120071 : iBean et al.ll2010D . Most of 
them are less dense than they would be if they were com¬ 
posed of silicate. Characterization of the relatively thick 
atmospheres of such exoplanets has been paid attention 
to as a first stepping stone towards an observational sci¬ 
ence as to habitable planets. 

A next stepping stone would be detection and char¬ 
acterization of thin atmospheres of close-in low-mass 
exoplanets, the observed mass-radius relationships for 
which are co nsistent with theor e tical ones for rocky 
planets (e.g., Fortnev et al.1120071 : IValencia et all 120071 : 
iGrasset et al'l l200iJh In particular, of special interest 
in this study are rocky super-Earths with masses of < 
10 M® and radii of < 2 i?® that are orbiting so close 
to their host stars that their surface temperatures are 
over the vaporization temperature of rock, because of 
intense stellar irradiation. Such exoplanets are re¬ 
ferred to as hot rocky super-Earths (HRSEs, hereafter) in 
this study. Examples o f the transi ting HRSEs detected 
so far are CoRoT-7 b ( jLeger et al l I2009f) , Kepler-10 b 
(|Batalha et all[201 111 . Kepler-78 b (Howard et al.l[20l3l) . 
and 55 Cnc e (TWinn et al.ll201lD . They are likely to have 
atmospheres with chemical species of rocky origin. 

Based on their chem ical equilibrium calculations, 
iSchaefer fe Feglevl (120091) presented the possible compo¬ 
sitions of the “silicate atmosphere” that is in gas-melt 
equilibrium with the molten rocky surface with no highly 
volatile elements such as H, C, N, S, and Cl (i.e., volatile- 
free magma ocean). They considered temperatures, T, 
of 1500-3000 K and gravity of 36.2 ms~ 2 , which cor- 
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responded to the planetary properties for CoRoT-7 b, 
and assumed several compositions for the volatile-free 
magma including the Earth’s continental and oceanic 
crusts, the bulk silicate Earth (BSE), and the bulk sil¬ 
icate Moon. Their calculations demonstrated that the 
atmosphere contained Na, O, O 2 , and SiO as the major 
constituents. Also, they investigated changes in atmo¬ 
spheric composition due to partial loss of the atmosphere, 
and found that Na was selectively lost from the atmo¬ 
sphere, while Mg became more abundant with fractional 
loss._ 

iMiguel et al.1 (I2011D made a similar but more exten¬ 
sive investigation of the composition of atmospheres of 
volatile-free HRSEs. They considered wider ranges of 
temperature (1000 to 3500 K), planetary mass (1 to 
10 M®), and planetary radius (l to 2.5 i?®). Also they 
considered komatiite in addition to BSE for the magma 
material. The resultant atm ospheric compositions ar e 
basically similar to those from ISchaefer fe Feglevl (I2009D : 
The major constituents are Na, O, O 2 , and SiO, although 
Fe and Mg are more abundant than SiO at T < 2000 K. 
They defined five types of HRSE atmosphere, depending 
on the abundance of the six species, and then classified 
Kepler planet candidates with radii of < 2.5 i?®, accord¬ 
ing to type. 

Since the presence of such atmospheres is specific to 
rocky planets, observational detection of those atmo¬ 
spheres could add a piece of convincing evidence for rocky 
planets. Furthermore, identifying the atmospheric con¬ 
stituents could give constraints on the bulk composition 
and formation process of HRSEs. However, the proper¬ 
ties and detectability of the atmospheres of the volatile- 
free HRSEs have not been examined yet. Thus, the aim 
of this study is to identify the radiative properties of 
the atmospheres on top of the global magma oceans of 
volatile-free HRSEs, and then assess their detectability. 
(In this study, we call such an atmosphere a mineral at¬ 
mosphere.) 

To this end, we model the thermal structure of the 
mineral atmospheres of HRSEs. The composition of the 
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atmosphere is obtained numerically in section [2] Also, 
we compile the absorption line data for the major gas 
species, which are published in the literature or online, 
and then calculate the absorption opacities in section [3] 
Using the chemical compositions and absorption opaci¬ 
ties, we integrate the thermal structure of the irradiated 
atmosphere in radiative equilibrium in section [I] Then, 
we evaluate the detectability of the chemical species in 
the mineral atmosphere through the emission spectra ob¬ 
tained during secondary eclipse in section [5] Finally, we 
discuss several effects to be examined in section [6] and 
then conclude this study in section [7] 

2. CHEMICAL COMPOSITION 
2.1. Element abundance and gas pressure 

To calculate the element abundance and pressure of 
gas in gas/melt chemical equilibrium with silicate melt 
for a given te mperature, we use the chemical equilibrium 
code MELTS (iGhiorso fc Sack! ri995t lAsimov fc Ghiorsol 
I1998D . which performs the free-energy minimization cal¬ 
culations. Figures [Q and [2] show the calculated total 
pressures and molar fractions of gas species, respectively, 
over a temperature range between 1500 and 3000 K 
for four magma compositions such as (a) BSE, (b) the 
Earth’s mid-ocean ridge basalts (MORB), (c) the Earth’s 
bulk crust, and (d) the Earth’s upper crust. The de¬ 
tail of the magma compositions is shown in Table [l] 
In this study, we co nsider the above f o ur co mpositions 
of magma, following iSchaefer fo Fcglevl (12009I ). although 
the magma composition of exoplanets is uncertain. The 
actual magma composition depends not only on the el¬ 
ement abundances but on the thermal history of the 
planet; exploring the details is beyond the scope of this 
study. Note that the choice among the four magma com¬ 
positions has little impact on the emission spectra in the 
infra-red wavelength region. Detectable difference may 
arise in the visible emission spectra, which is discussed 
in section f6. 61 

The overall trend shown in Fig s. [I] and [2] i s similar with 
that fo und bv lSchacfer fc Feglcvl (I2009T ) and lMiguel et al.l 
Umi), who used the MAGM A chemical equilibrium 
code dFeglev fc Cameronl 119871) . instead of MELTS. For 
example, Na is the most abundant species in most of the 
temperature range, while SiO increases with temperature 
and becomes the most abundant one for T > 2800 K for 
BSE and MORB, T > 2600 K for the bulk crust, and 
T > 2500 K for the upper crust. Also, the total vapor 
pressure is as small as ~ 1CU 7 bar at T = 1500 K and 
~ 0.1 bar at T = 3000 K, for example (see Fig.[l]). 

A noticeable difference is seen in the abundance of 
potassium gas be t ween our calculations and tho se by 
ISchaefer fc Fegletl (120091) and iMiguel et al.1 (1201 111 . Our 
calculations yield much higher partial pressure of potas¬ 
sium: The ratio of the potassium to sodium partial 
pressures that we have obtained is larger by a fac- 
tor of more than 100, relative to that calculated by 
ISchaefer fc Fegletl (.2009). This difference is due mainly 
to that in the adopted thermodynamic model for sil¬ 
icate melt. The melt model of MAGMA used in 
ISchaefer fe Feglevl (120091) is ideal mixing of many sim¬ 
ple and complex fictitious oxide species, which is cali¬ 
brated for silicate melts with high-K 2 0 contents, whereas 
MELTS used in this study adopts a symmetric regular 


solution model with minimum melt components, which 
is calibrated by experimental results on natural systems 
with I 0 W-K 2 O abundances. MAGMA would appropri¬ 
ately estimate potassium vapor pressure fo r high -K?Q 
melt as shown in Fig. 5 of lSchaefer fc Feglevl (120041) . but 
is not good for I 0 W-K 2 O melt. By contrast, MELTS gives 
appropriate results for I 0 W-K 2 O melt, which is of interest 
in this study (see Table [TJ. The impact on the emission 
spectrum from the mineral atmosphere is discussed in 
section 16.61 


Table 1 

Bulk compositions considered in this study 


Oxide(wt%) 

BSE“ 

MORB” 

Bulk crust 6 

Upper crust 6 

sio 2 

45.1 

49.6 

60.6 

66.6 

MgO 

37.9 

9.75 

4.70 

2.50 

FeO 

8.06 

8.06 

6.71 

5.04 

ai 2 o 3 

4.46 

16.8 

15.9 

15.4 

CaO 

3.55 

12.5 

6.41 

3.59 

Na 2 0 

0.36 

2.18 

3.07 

3.27 

Cr 2 0 3 

0.38 

0.07 

0.0004 

0.0003 

Ti0 2 

0.20 

0.9 

0.72 

0.64 

k 2 o 

0.03 

0.07 

1.81 

2.80 

p 2 o 5 

0.02 

0.10 

0.13 

0.15 


Notes. FlMcDonough fe Sun! ll 1 9951 ) . PlIH iitinick fc Gad 120031 1 



Figure 1. Total pressures of gas in gas/melt equilibrium with 
magma of assumed composition are shown as functions of temper¬ 
ature. The assumed bulk compositions are the same as those of 
the bulk silicate Earth (BSE) (red), the Earth’s mid-ocean ridge 
basalts (MORB) (green), the Earth’s bulk crust (blue), and the 
Earth’s upper crust (magenta) (see Table [T] for the details). 


2.2. Gas phase equilibrium 

To calculate the equilibrium abundances of gas species 
i n the m in eral atmos p here, we use the NASA CEA code 
(I Gordon fe McBridel I1996D . which performs the free- 
energy minimization calculations. In practice, we cal¬ 
culate the abundances for some selected temperatures 
from 2000 K to 4500 K and pressures from 10~ 8 bar to 
10 _1 bar, which are of interest in this study, consider¬ 
ing the gas phase equilibrium of a system composed of 
only the major elements such as O, Na, K, Si and Fe. 
We assume that the element abundances are vertically 
constant; their values are calculated in section l2Jl 
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(a) BSE 


(b) MORB 



(c) Bulk crust 



(d) Upper crust 



Figure 2. Composition of gas in chemical equilibrium with molten silicate. Molar fractions of the gas species are shown as functions of 
temperature. The assumed bulk compositions are the same as those of (a) the bulk silicate Earth (BSE), (b) the Earth’s mid-ocean ridge 
basalts (MORB), (c) the Earth’s bulk crust, and (d) the Earth’s upper crust (see Tablc[l]for the details). 


In our atmosphere model, the molar fractions of Na, K 
and Fe gas vary little, compared to those of O, 0 2 , Si and 
SiO. 0 2 and SiO tend to decrease and O and Si increase 
with altitude in low pressure regions in the atmosphere 
(see section a, because of thermal dissociations of 0 2 
and SiO. 

Although we ignore photo-chemical reactions in this 
study, the photo-dissociation of SiO may be important 
for the detectability of the atmosphere. We discuss the 
possibility of the photo-dissociation of SiO and the im¬ 
pact on the atmospheric structure and detectability in 
section 16.31 


3. ABSORPTION OPACITY 

Among the chemical species shown in Figure [2] we 
consider line absorption by seven major species that in¬ 
clude Na, K, Fe, O, 0 2 , Si and SiO. We neglect Rayleigh 
scattering, because the Rayleigh scattering opacities are 
mostly smaller than the absorption opacity of SiO at 
wav eleng ths of interest in this study, as discussed in sec¬ 
tion [63] Note that whil e TiO is known to be a stron g 
absorber for visible ray ([Allard et al.11200(1 iPlezI 119981 ). 
its absorption opacity is always smaller than Na absorp¬ 
tion opacity in the visible and IR wavelength regions, be¬ 
cause the molar fraction of TiO is negligibly small (e.g., 
- 10" 12 at T = 1500 K and ~ 10“ 6 at T = 3000 K), com¬ 
pared to that of Na (> 0.1) in the atmosphere considered 


here. 

Absorption opacities are calculated in the fol¬ 
lowing way, except for 0 2 , for which we ap¬ 
ply vibrat ional-rotational trans itions in the HITRAN 
database ([Rothman et al.l I2013D . In local thermody¬ 
namic equilibrium, the absorption cross section of species 
A for an energy transitio n from the zth level to the jth 
level, is given as ([Piskunov fc Kimkai 12001 ) 





7re 2 exp(— Ei/kT) 


m e c 
1 — exp 


Q(T) 

'Ej-EA 


H(u, a), 


( 1 ) 


where v is the frequency, gi is the statistical weight of 
the zth energy level, / ;J - is the oscillator strength for 
a transition from the ith (lower) to jth (higher) level, 
Ei and Ej are excitation energies of the ith and jth 
energy levels, respectively, and Q(T) is the partition 
function of the species at temperature T. e is the ele¬ 
mentary charge (= 4.803 x 10 -10 esu), m e is the elec¬ 
tron mass (= 9.1094 x 10” 28 g), c is the light velocity 
(= 2.9979246 x 10 10 cms -1 ), and k is the Boltzmann 
constant (= 1.3806488 x 10~ 1 8 ergK" 1 ). W e take the 
values of q,;, / , an d Ei from iKurucd ([1992D (see also 
iBarton et al.l ( 20131 ) for SiO). Also we take the values 
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Wavelength [nm] Wavelength [jaiti] 


Figure 3. Absorption cross sections of the seven major gas species, SiO (red), O 2 (green), Na (blue), K (magenta), Fe (grey), O (cyan), 
and Si(yellow) as functions of wavelength. Each cross section is multiplied by each molar fraction. The magma composition is assumed 
to be BSE (see Table [Q. (a): The temperature is 2000 K and the pressure is 1 x 10 —5 bar. The calculated molar fractions of Na, K, Fe, 
Si, SiO, O, and O 2 are 0.721, 0.0407, 0.001, 1 x 10 —6 , 0.00731, 0.0233, and 0.189, respectively, (b): The temperature is 3000 K and the 
pressure is 1 x 10 —2 bar. The calculated molar fractions of Na, K, Fe, Si, SiO, O, and O 2 are 0.204, 0.0142, 0.01, 1 x 10 —3 , 0.343, 0.0991, 
and 0.211, respectively. 


of Q(T) for Na, K, Fe, S i and O from llrwliil (|1981h and 
for SiO from iRossi et all (Il985li . 

The f unctio n H(u,a) is the l ine p rofile function given 
by (e.g.. [Rvbicki fe Lightmanlll986l) 


H(u , a) 


a ) 

s/nAv D ' 


( 2 ) 


where a is the line width parameter, V(v, a) is the Voigt 
function, and A ud is the Doppler width given by 


2kT 


Av d = ^\—, (3) 

c V m A 

where tua is the mass of the absorber and is the 
central frequency of the spectral line. The line width 
parameter a is calculated with damping constants for 
natural broadening and van der Waals broadening as 
dPiskunov fe Kupkall200lD 

a = (4) 

_47rA pn 

According to lUnsoldl (|1955l) and ICravl (119761) . T r and T w 
are given by 

T r = 0.222 x 10 14 A" 2 , (5) 


T\y = 17 X 


^0.3 x 1CT 30 


' (* + l) 2 
.{i-Xi? 


x v 3/5 N, 


(1 + l ) 2 
C i-Xi ) 2 

( 6 ) 


where A is wavelength in nm, N is the total number 
density of the gas species in cnr 3 , / is the ionization 
potential in eV, \i is the excitation potential at the ith 
level in eV, and l is the electric charge, v is the relative 
thermal velocity defined by 


v = v 


A 



(7) 


where v A is the thermal velocity of the absorber and [i A 
and p are the molecular weight of the absorber and the 
mean molecular weight, respectively. As for the SiO ab¬ 
sorption opacity, we ass ume Tw = 10 -7 x N , following 
iKurucz fe Avrettl (11981D . To calculate the line profile 


H (i/ : a) numerically, we us e a polynom ial expansion of 
the Voigt function given bv iKimty] (11997D . We have con¬ 
firmed that this approximation yields relative errors of 
the order of 1(R 4 , which is small enough for our discus¬ 
sion in this study. 

Figure [3] shows the cross section of each species multi¬ 
plied by each molar fraction as a function of wavelength 
for T = 2000 K(a) and 3000 K(b), respectively. The as¬ 
sumed pressures are 1.0 x 10~ 5 bar for T = 2000 K and 
0.01 bar for T = 3000 K. The magma composition is as¬ 
sumed to be BSE. The calculated molar fractions of Na, 
K, Fe, Si, SiO, O, and O 2 are, respectively, iN a = 0.721, 
x K = 0.0407, x Fe = 0.001, z S i = 10“ 6 , z S iO = 0.00731, 
xo = 0.0233, and xq 2 = 0.189 for T = 2000 K(a) and 
iNa = 0.204, x K = 0.0142, x Fe = 0.01, x Si = l x 10" 3 , 
^SiO = 0.343, xo = 0.0991, and xo 2 = 0.211 for T = 
3000 K(b). 

In Fig. [3K, the most abundant species Na (blue) 
presents a few hundred of strong absorption lines between 
0.253 and 10.0 pm, including the D line at 0.5893 pm. 
Also, several hundred strong lines of K (magenta) that 
are distinguishable from the Na lines are also found be¬ 
tween 0.286 and 10.0 pm, including potassium doublet 
lines at 0.766 and 0.770 pm. While SiO (red) is less 
abundant by a factor of about 100 than Na, millions of 
its electron-transition absorption lines dominate in the 
FUV wavelength region (0.14 - 0.3 pm) and its absorp¬ 
tion due to rotation-vibration transition dominates in 
several parts of the IR wavelength region, except narrow 
line peaks of Na, K and Fe. The O electron-transition 
absorption (cyan) dominates between 0.10 and 0.13 pm. 
The several thousands lines of Fe (grey) dominates be¬ 
tween 0.20 and 0.54 pm. Although O 2 (green) is the 
second most abundant species, its absorption is obscured 
by the other absorptions except at a narrow range around 
30 pm. Also, the Si electron-transition absorption (yel¬ 
low) is obscured by the SiO absorption in both Figs. [3ji 
and [3}o and presents many prominent absorption lines 
between 0.1 and 0.3 pm. 

In the higher-T case shown in Fig. [3^, the SiO fea¬ 
ture is more prominent. The electron-transition absorp¬ 
tion in the FUV region is larger by about two orders 
of magnitude than that in Fig. [5|i, because the molar 
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fraction of SiO is larger by two orders of magnitude (see 
Fig. © and comparable with Na at 3000 K. The most re¬ 
markable difference in absorption feature between Fig.[3[i 
and Fig. [3jo is that the SiO rotation-vibration absorption 
dominates at all the wavelengths longer than about 3 /mi, 
except several hundreds of narrow lines of Na and K. The 
rises in the peak values of SiO are due to the increase in 
abundance, while the enhancement of pressure broaden¬ 
ing also contributes to the rises in absorption between 
those peaks. 

Other minor species might contribute to the opacity. 
As for Mg and P present in Fig. [21 we have confirmed 
that the absorption opacities of these species are negligi¬ 
bly small relative to the SiO absorption opacity in most 
of the wavelength region of interest, except for narrow 
ranges between 0.1 and 0.15 /im. Absorption line data 
of the diatomic molecules other than SiO in Fig. [2] are 
unavailable, to our knowledge. Thus, we have ignored 
their contribution. 

Finally, for integrating the thermal structure of the at¬ 
mosphere in the next section, we use the total absorption 
opacity, calculated as 

pK v = N^2x A cr^(i, j), (8) 

A 

where p is the mass density. 

4. TEMPERATURE PROFILE 

In this study, we simulate the one-dimensional plane- 
parallel thermal structure of the atmosphere that is in the 
radiative, hydrostatic, and chemical equilibrium. To do 
so, we integrate the so-called two-stream equations with 
the assumption of quasi-isotro pic radiation, ado pting the 
S —Eddington approximation dToon et al.l[1989 ): 

+ i K-wn m 

^- = -\F+- 7 -F-+2nB v (T), ( 10 ) 

where F+ and F~ are the upward and downward radia¬ 
tive fluxes, respectively, B„{T) is the Planck function, 
and r„ is the optical depth defined by 

dr v n i \ 

— = -K v p, (11) 

az 

where 2 is the altitude. In equations © and m , we 
have assumed that the single-scattering albedo is 0, be¬ 
cause absorption dominates scattering in most of the 
wavelength region of interest, as described above. The 
boundary conditions at the top and bottom of the at¬ 
mosphere are respectively F~ = 0 and F+ = ttB v (T^), 
where T g is the ground temperature. To solve equa¬ 
tions © and cnD, we have used the w idely-used al¬ 
go rithm devel o ped b y iToon et all (119890 . According 
to IToon et ahl (119890 , the error in the radiative fluxes 
caused by the two-stream approximation was typically 
less than 10%, which is small relative to uncertainties in 
physical properties of exoplanets. 

The net radiative flux integrated over all frequencies 
is constant throughout the atmosphere in radiative equi¬ 
librium. For an irradiated atmosphere, the net radiative 
flux, F net)I/ , is written as 

Fnet,v(r v ) = F+(r„) - F~(t v ) - /r*F„*e- T ^*, (12) 


where F* is the incident stellar flux and /r* is the cosine 
of the zenith angle of the incident stellar radiation. The 
radiative equilibrium condition is given as 

pOO 

/ Faet,u dlS = Fq, (13) 

JO 

where Fq is the constant flux. An HRSE is subject to 


tidal heating. Its strength de pends on orbital a nd in¬ 
ternal properties of the HRSE ([Barnes et al.lf2010h . The 
detailed consideration of the tidal heating effect is be¬ 
yond the scope of this study (A brief discussion is made 
in section 16.411 . In this study, we assume Fq = 10 W/m 2 
for numerical stability. For Fq < 10 W/m 2 , we have 
confirmed that choice of Fq have little influence on the 
thermal structure of the atmosphere. 

We consider a pressure range from 1 x 10~ 8 bar to the 
pressure at ground level, P g . Note that P g is determined 
a posteriori, because P g is the vapor pressure, which de¬ 
pends on the ground temperature, T g . Numerically, the 
atmosphere is vertically divided into 50 layers. The lay¬ 
ers are prepared so that the size of each layer logarith¬ 
mically increases with pressure. In addition, to integrate 
equations © and (flOl) . we consider 100 spectral inter¬ 
vals in a frequency range between lxlO 13 and 2xl0 15 
s -1 such that the size of each interval logarithmically in¬ 
creases with frequency. We calculate the harmonic mean 
opacities using cross sections obtained in section [3] The 
calculation method of the mean opacities is described in 
Appendix [A] 

To find the thermal structure of the atmosphere in ra¬ 
diative equilibrium and also in gas-melt chemical equi¬ 
librium with the underlying magma ocean, we need 
three iterative procedures. Firstly, to find the radiative- 
equilibrium solution for given opacities and molar frac¬ 
tions of atmospheric gases, we integrate the equation of 
energy conservation, 


dT n _ dF n 
dt ~ dP n ’ 


(14) 


where t is a normalized time, T n , P n and F n are the 
temperature, pressure, and net radiative flux in the nth 
layer, respectively. The integration is continued until 
| F n — Fol /Fq becomes smaller than 1 % in all the layers. 

Secondly, once we find the radiative-equilibrium solu¬ 
tion, we calculate the opacities that are consistent with 
temperatures obtained above. In practice, to save mem¬ 
ory and epu time, we have prepared a numerical table in 
which the harmonic mean opacities (see Appendix [Aj in 
100 spectral intervals are given as functions of T, log P, 
and, the mean molecular weight, Ji\ The grids are pre¬ 
pared for T = 2000, 2500, 3000, 3500, 4000, and 4500 K, 
log(P/bar) = —1, —2, —3, —4, and —5, and JI = 25, 
30, and 35. In the pressure range from 1 x 10 -8 bar to 
1 x 10 -5 bar, we have ignored the pressure dependence 
of the opacity because r r > T w in the range. Note that 
we have checked that the calculated flux differs by at 
most 5 %, even if a smaller table with T = 2000, 3000, 
and 4000 K, log(P/bar) = —1, —3, and —5, and JL = 25 
and 35 is used. The convergence condition is that the 
temperature in each layer does not vary within less than 
1 %. 

Finally, we have to find the solution in which the at¬ 
mosphere is in equilibrium with the underlying magma 
ocean at its bottom. The ground pressure P g and mo¬ 
lar fractions xa are functions of the ground temperature 
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Figure 4. Temperature-pressure profile in the mineral atmo¬ 
sphere on top of the BSE magma of a super-Earth with gravity of 
25 m/s 2 for five choices of the substellar-point equilibrium temper¬ 
ature, T eq , of 1800 K (black), 2000 K (cy an), 2300 K (blue), 2500 K 
(green), and 3000 K (red) (see eq. 1161 ; A p = 0). The host star 
is assumed to be a Sun-like star with radius of 1 Rq and effective 
temperature of 6000 K and emit the blackbody radiation of 6000 K. 
The solid lines show the day-side averaged profiles (i.e ., the cosine 
of the stellar-light zenith angle //* = 1/2; see eq. 1121 '). The open 
circles show the temperatures at the bottom of the atmosphere, 
7b, and the filled circles show the temperatures at the ground, T g . 
The orange dotted line represents the total vapor pressure for the 
BSE composition, which corresponds to the ground. 

T g (see section (3). From equations and (fl3l) , the 
magma being assumed to be a blackbody, T g is given as 

aT* = F 0 + J^ {f~ v {P s ) + dv, (15) 

where a is the Stefan-Boltzmann constant (= 5.67 x 
lCT 8 W/m 2 K 4 ) and r v>g is the total optical depth of the 
atmosphere. We calculate T g from Eq. m and then cal¬ 
culate P g and Xa • The gas chemical equilibrium compo¬ 
sition in each layer, is obtained by interpolating lin¬ 
early the tabulated data made by the NASA CEA code. 
The grids are prepared for T = 2000, 2500, 3000, 3500, 
4000, and 4500 K, log(P/bar) = -1, -2, -3, -4, -5, -6, 
—7 and —8, and T g being every one kelvin from 1500K 
to 3000K. The convergence condition is that T g and the 
temperature in each layer do not vary within less than 
1 %. 

Figure H] shows the calculated temperature profiles for 
different substellar-point equilibrium temperatures de¬ 
fined by 

T 4 q =(l-A P )|§T 4 , (16) 

where 1?* and T* are respectively the radius and temper¬ 
ature of the host star, A p is the planetary albedo, and D 
is the orbital distance of the planet. Five values of T eq 
are chosen: T eq = 1800 K (black), 2000 K (cyan), 2300 K 
(blue), 2500 K (green), and 3000 K (red). The incident 
stellar flux is calculated as F* = T 4 q /T 4 x nB u (T :t ). The 
solid lines represent the temperature profiles at p,» = 1/2 
(i.e., day-side average). Also, the open and filled cir¬ 
cles show the temperatures at the bottom of the atmo¬ 
sphere, Tb, and at the ground, T g , respectively. Figure [5] 
shows the calculated molar fraction profiles of seven ma¬ 
jor species for the atmosphere structure shown in Fig. Q] 
Here we have assumed g = 25 m/s 2 and A p = 0. The host 
star is assumed to emit blackbody radiation of 6000 K. 


(a) T eq = 1800K 



(b) Teq = 2000K 



(c) Teq = 2300K 



(d) Teq = 2500K 



(e) Teq = 3000K 



Figure 5. The distribution of the seven major gas species in 
chemical equilibrium in the mineral atmosphere. Molar fractions 
of SiO (red), O 2 (green), Na (blue), K (magenta), Fe (grey), O 
(cyan) and Si (yellow) are shown as functions of pressure for five 
choices of the substellar-point equilibrium temperature; (a) T eq 
= 1800K, (b) 2000K, (c) 2300K, (d) 2500K, and (e) 3000K. The 
corresponding temperature-pressure structure of the atmosphere is 
shown in Figure l4l 

The magma composition is assumed to be BSE. 

The most remarkable feature of the temperature pro¬ 
files in Fig. |4] is the thermal inversion structure in the 
cases of T eq > 2300 K. As found in section [3l the ab¬ 
sorption is more significant in the UV and visible wave¬ 
length regions (0.1 to ~1 /im) than in the IR region (~1 
to 4 /im), mainly because of the strong UV-ray absorp¬ 
tion by SiO and visible-light absorption by the doublet 
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lines of Na and K. (Note that most of the absorption 
lines of Fe, Na and K are too narrow to contribute to 
the mean opacity.) This means that the absorption 
of the incident stellar radiation is stronger compared to 
the intrinsic planetary radiation, which results in ther¬ 
mal inversion for P > 10~ 5 bar. For P < lCP 5 bar, the 
temperature decreases with decreasing pressure, because 
the dissociation of SiO occurs as shown in Figure 0J so 
that the UV-ray absorption by SiO becomes weak. 

Not only the thermal inversion region extends but T g 
increases, as T eq increases. In the case of T eq = 3000 K, 
T g , is significantly smaller than /x* T eq . In contrast, 
T g ~ /x* T eq for low T eq , because the atmosphere is so 
optically thin that the ground is directly heated by stellar 
irradiation. In the cases of T eq < 2000 K, the atmosphere 
is thus isothermal. Note that jumps in temperature be¬ 
tween Tb and T g are found in Figure 0] This is because 
those atmospheres are optical thin for the incident stellar 
flux and any heat transfer process except radiative trans¬ 
fer is not taken into account. Thus, in reality, the at¬ 
mosphere is not exactly isothermal, because there should 
be a conductive region with steep temperature gradient 
near the ground. Such structure, however, has little 
impact on the emission spectra that we use when dis¬ 
cussing the detectability of HSREs in secondary eclipse 
(see section 0. 


5. DETECTABILITY VIA SECONDARY ECLIPSE 
OBSERVATION 


5.1. Secondary eclipse depth 

Here we assess the detectability of the mineral atmo¬ 
spheres of HRSEs that we have modeled via the sec¬ 
ondary eclipse observation. This observation enables us 
to infer the vertical temperature profile and composition 
of an exoplanetary atmosphere by obtaining planetary 
emission spectrum. In particular, the thermal inversion 
structure found in section 0 ] is expected to be detected. 

The secondary eclipse observation measures the ratio 
of the planetary luminosity to the sum of the stellar and 
planetary luminosities. The ratio is often called the sec¬ 
ondary eclipse depth, ex , which is given by 


/ R P \ F Py \ + ApF*^ 
\R*) T*^ 


(17) 


where F. t \ is the emergent stellar flux and F Pt \ is the 
emergent planetary flux calculated as 

nOO 

F p ,\ = 7r / B\(T(r\))e~ Tx dT\. (18) 

Jo 

Below, as in section 01 we assume g = 25 m/s 2 , A p = 0, 
and blackbody stellar radiation of 6000 K. The magma 
composition is assumed to be BSE. Also, we suppose an 
HRSE with radii of 2 orbiting a solar analog, namely 
Rp/R* = 0.02. To integrate equations (ITTl) and ([T5|h 
we consider 30000 spectral intervals in the wavelength 
range between 0.1 and 100 /xm; the size of each interval 
logarithmically increases with wavelength. Then, as in 
section 0J we calculate the arithmetic mean opacities, 
interpolating linearly the tabulated data. We consider 
only the case of /x* = 1/2, because we are interested in 
planetary disk averaged emission (i.e., day-side average). 

Figure [6] plots the secondary eclipse depth e\ (and also 
F p ,x/F^x) for T eq = 1800 K (black), 2000 K (cyan), 


o 



Je a =3000K - 
TS=2500K 
TS=2300K 
Tpn=2000K 

Teq = 1800K 


1 10 
Wavelength [pm] 



Wavelength [pm] 


Figure 6. Predicted dayside-averaged emission spectra of a hot 
rocky super-Earth of 2 Rq in secondary eclipse that has a mineral 
atmosphere in equilibrium with an u nde rlying BSE magma ocean. 
The secondary eclipse depth (see eq. EH) is shown as a function of 
wavelength in the range of (a) 0.1-100 pm and (b) 0.3-1 pm. The 
inset in the Fig. (a) is an enlarged view of the spectra in 3-5 pm. 
Five equilibrium temperatures are chosen: T eq = 1800 K (black), 
2000 K (cyan), 2300 K (blue), 2500 K (green), and 3000 K (red). 
The corresponding atmospheric structures are shown by dashed 
lines in Fig. [4] The host star is assumed to be a Sun-like star with 
radius of 1 R@ (i.e., the planet ary/stellar radius ratio being 0.02) 
and emit the blackbody radiation of 6000 K. 


2300 K (blue), 2500 K (green), and 3000 K (red): Pan¬ 
els (a) and (b) show the spectra in the wavelength re¬ 
gion of 0.1-100 /xm and 0.3-1 /xm, respectively. The 
secondary-eclipse-depth spectra for T eq = 1800 K and 
2000 K exhibit only some narrow line features of Na and 
K. This is mainly because the atmospheres are isothermal 
(see Fig. 01 ) except for the temperature jump above the 
ground. The line features are, thus, due to the emission 
from parts of the atmosphere with temperatures higher 
than Tg, while the continuous spectra are the black body 
spectra of T g . 

By contrast, the spectra for T eq = 2300 K and 2500 K 
present noticeable features of SiO around 10 and 100 /xm, 
in addition to hundreds of narrow features of Na and K. 
For T eq = 3000 K, a spectral feature by SiO around 4 /xm 
appears in addition to those around 10 and 100 /xm, as 
shown in the subfigure inside Fig[6}i. The features of 
SiO become stronger as T eq increases, because the SiO 






































































































































abundance is sensitive to temperature and also the ther¬ 
mal inversion is more significant at higher temperature. 
Additionally, one can also identify features of the sodium 
doublet at 0.589 fim and the potassium doublet at 0.77 
/jm and narrow line features of Fe, Na and K in the vis¬ 
ible band, as shown in Panel (b). 


are assumed to be 100 in 0.1-1 /im and 10 in 1-100 fim. 
affect the value of (S/N)ax- As shown in Figure [Tji, 
the spectral features of SiO are sufficiently detectable at 
4, 10 and 100 fim for T eq = 3000 K. The SiO feature is 
also sufficiently detectable at 10 fin i for T eq = 2500 K 
(Fig. [TJd), while the feature is marginally undetectable 


5.2. Detectability 

We evaluate the detectability of the mineral at¬ 
mosphere of an HRSE for near-future space-based 
secondary-eclipse observation. The photometric accu¬ 
racy of the transit/eclipse observation from space de¬ 
pends on the instrument measurement noise, the intrinsic 
stellar va riabi lit y, and the sho t noise from photon statis¬ 
tics (e.g. iClanton et al.l 120121) . In this section, we quan¬ 
tify the photometric accuracy required for the line detec¬ 
tion and discuss the feasibility, assuming the ideal case 
that the total noise is dominated by the shot noise (i.e. 
the photon-noise limit). To quantify the detectability for 
the photon-noise limit, we introduce the signal-to-noise 
ratio of the detection of the secondary eclipse as 

<S/iVW = J=, (19) 


_ , , / Atobs -R* 7tRaa(F*) /on', 

—i —• <20) 

where AA is the band width used for the detection of the 
secondary eclipse, IV* is the photon counts of the star, A 
is the photon collecting area, t Q b s is the exposure time, d 
is the distance of the system, 91 is the resolving power of 
the spectra, ca,aa is the depth of the secondary eclipse 
within AA, N p is the photon count of the planet, Raa(T*) 
is the blackbody radiation integrated over a wavelength 
band (i.e., J AX B\(T*)dA/A\), and h is the Planck con¬ 
stant. 

Figurc[7]plots the theoretical (solid line), mock (cross), 
and blackbody (dotted lines) secondary-eclipse-depth 
spectra for T eq = 3000 K (a), 2500 K (b), and 2300 K (c). 
When calculating the mock spectra, we us e the false nor¬ 
mal random number generation circuit dWilliam ct al.l 
11996 1. and include the statistical noise, assuming the 
normal distribution with the standard deviation a = 
VK- We assume A = (5/2 m) 2 as a fiducial value, 
roughly correspon ding to that o f James Webb Space 
Telescope (JWST; iClampinl 120081) . and that 91 = 100 
in 0.3-1 fin i and 10 in 1-100 fim, i 0 b s = 10 hours, and 
d = 100 pc. The blackbody spectra are calculated from 
equation na with the assumption F p \ = B\{T). 

As demonstrated in Fig. [71 the spectral features of Na, 
K and SiO are represented by the blackbody spectra with 
different temperatures (i.e., brightness temperature Tb r ). 
For example, in the case of T eq = 3000 K (see Fig.jT^i), Tbr 
~ 3200 K at 10/rm, and 2700 K at 4/xm due to SiO. Ad¬ 
ditionally, although not shown, Tb r — 3500 K at 0.6 fim 
and at 0.8^m (which are due to Na and K, respectively), 
and 3100 K at 100 fim (SiO). At the other wavelengths, 
Tbr — T g (e.g., Tbr = 2350 K for T eq = 3000 K), because 
the atmosphere is optically thin at those wavelengths. 

Comparing the mock spectra with the blackbody spec¬ 
tra, one can quantify the detectability of the spectral fea¬ 
tures of Na, K and SiO. The error bars become larger in 
both shorter and longer wavelength regions, because the 
number of incident photons (i.e., B aa(T*)) is small in 
such wavelength regions. Also, the values of 91, which 
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Figure 7. The secondary eclipse depths theoretically predicted 
and the mock observational spectrum compared to the blackbody 
spectra (BB). The secondary eclipse depths are shown as functions 
of wavelength for T eq = (a) 3000 K, (b) 2500 K, and (c) 2300 K. 
The black dotted lines show the secondary eclipse depths of the 
blackbody, temperatures of which are indicated by ”BB(T)”. The 
bars show the mock spectrum calculated based on the photon- 
noise limit with the assumptions that R * = 1 Rq, R p = 0.02 Rq, 
A = 7r(5m/2) 2 , £ 0 b s = 10 hours, = 100 in 0.3-1 fim and 10 
in 1-100 fim, and d = 100 pc (see the text for the details). The 
corresponding atmospheric structures are shown by dashed lines in 
Fig- GO 
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Figure 8. (a) The secondary eclip se d epth relative to that of 

blackbody with T g ( rj\ aa; see eq. 1211 ) and (b) the minimum 
signal-to-noise ratio required to distinguish the emission spectral 
feat ures from the blackbody emission with T g (|(5'/-/V)AAlmin; see 
eq. 1221 ). 51 = 100 in 0.3-1 pm and 10 in 1-100 pm. Three equi¬ 
librium temperatures are chosen: T eq = 3000 K (red), 2500 K 
(green), 2300 K (blue). The corresponding atmospheric structures 
are shown by dashed lines in Fig. [4] The host star is assumed to be 
a Sun-like star with radius of 1 Rq (i.e., the planetary/stellar ra¬ 
dius ratio being 0.02) and emit the blackbody radiation of 6000 K. 

for T eq = 2300 K (Fig. Eh). 

The signal-to-noise ratio for detecting the line absorp¬ 
tion spectral feature of the secondary eclipse depth is 
generally smaller than that for the secondary eclipse de¬ 
tection. Thus, we define the secondary eclipse depth rel¬ 
ative to that with no atmospheric feature (i.e., blackbody 
with T g ), 7 ?a,aa as 


Pa, aa 


r \+A\ 


£A 


dX 

AA 


( 21 ) 


where £BB,A(T , g ) is the secondary eclipse depth by the 
blackbody radiation with T g . Using p a,aa ; we estimate 
the detectability of the line spectral features. From 
equations and m, we introduce a new quality, 
(S/N)l : aa, that is defined as 

(S/N) LAX =rj XAX (S/N) AX . (22) 

When ( S/N)l AX is sufficiently large, we can detect 
the line spectral feature of the secondary eclipse depth. 
Thus, if the minimum value of {S/N)l Ax required 
to detect a line spectral feature, |(iS'/-/V)i /i AA|min, is 3, 



Figure 9. The limiting distance d max (see eq. 1231 for the def- 
inition), interior to which an HRSE with the model atmosphere 
orbiting a Sun-like star is detectable with (S/N)/±\ by a given 
instrument, is shown as a function of wavelength for three equi¬ 
librium temperatures T eq = 3000 K (red), 2500 K (green), and 
2300 K (blue) with the assumption that R* = 1 Rq, R p = 0.02 Rq, 
A = 7r(5m/2) 2 , t Q b s = 10 hours, and 5t = 10 in 1-100 pm and 100 
in 0.3-1 pm. The corresponding atmospheric structure is shown in 
Fig. a The solid lines show the limiting distance d max to detect 
the line spectral features within secondary eclipse depth around 4, 
10, and 100 pm, assuming ( S/N)/\\ = 3/77 shown in Figure 03. 
The dotted lines show the limiting distance cf max to detect the sec¬ 
ondary eclipse depth itself, assuming (S/N)/±\= 3. The host star 
is assumed to be a Sun-like star with radius of 1 Rq (i.e., the 
planet ary/stellar radius ratio being 0.02) and emit the blackbody 
radiation of 6000 K. 


{S/N) aa > 3/?7a,aa is required to do so. 

Figure [HJi shows r) XAX for T eq = 3000 K (red), 2500 K 
(green), and 2300 K (blue) with SR = 100 in 0.3-1 pm 
and 10 in 1-100 pm. As shown in Figure [8^,, pa.aa takes 
large values at some wavelengths, because of the strong 
line features: r) XAX ~ 1 at 0.6 pm (Na) and ~ 0.35 at 
10 pm (SiO) for T eq = 3000 K (red), ija.aa is large at 
~ 4 pm because of the SiO feature peak only for T eq = 
3000 K (red), since the SiO feature is weaker than Na and 
K features for T eq = 2300 K and 2500 K. Also, Figure [8 ]d 
shows | {S / V)^ ; AA |min/ PA,Aa(— \ (S/N) AX |min) for T eq = 
3000 K (red), 2500 K (green), and 2300 K (blue), assum¬ 
ing | (5/V) l AX |min = 3 as a minimum fiducial value to 
detect the line spectral feature with SR = 100 in 0.3-1 pm 
and 10 in 1-100 pm. We obtain that |(,S/iV)AA|min < 100 
around 4, 10 and 100 pm; for example, |(S/lV)AA|min — 
60 at 10 pm for T eq = 2300 K (blue) and |(S/lV)AA|min — 
10 at 4 pm for T eq = 3000 K (red). Thus, observation of 
an HRSE in secondary eclipse can detect the components 
of the mineral atmosphere such as Na, K and SiO, pro¬ 
vided its signal to noise ratio is larger than |(S/iV)AA|min- 

Finally, we estimate the limiting distance d max interior 
to which an HRSE with the mineral atmosphere is de¬ 
tectable with a given |(S'/A r )AA|min, for the photon-noise 
limiting case. From equation (12011 . G? max is written as 

dm ax {|(S/A0AAUin} = . 

\{S/JX ) AA | min 


X 


At 0 bs tHaa(F*) 

Tr h 


R*. (23) 


Figure [9] shows d max as a function of wavelength for 
| (SyV) aa| min = 3 /pa.aa (he., the minimum signal-to- 
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noise ratio for the detection of line spectral features 
within secondary eclipse depth; shown by solid lines) or 3 
(i.e., the minimum signal-to-noise ratio for the detection 
of the secondary eclipse depth itself; shown by dotted 
lines); R* = 1 Rq, R p = 0.02 Rq, A = 7r(5m/2) 2 , t obs = 
10 hours, and 91 = 10 in 1-100 /im and 100 in 0.3-1 /im. 
In particular, the solid lines show the limiting distance 
dmax to detect the SiO line spectral features around 4, 
10, and 100 /an. Note that G? max around 4 /tm is not 
shown for T eq = 2300 K and 2500 K, because the Na, 
K and SiO features in those cases are indistinguishable 
from each other when integrating the secondary eclipse 
depth with respect to wavelength with 91 = 10. 

Table [2] summarizes d max for the SiO detection. The 
limiting distance d max increases with increasing T eq , be¬ 
cause an HRSE with higher T eq shows greater planetary 
emission and more remarkable line absorption spectral 
feature of the secondary eclipse depths. Thus, the min¬ 
eral atmospheres of more distant HRSEs would be also 
detectable if they have higher T eq . 


Table 2 

Limiting distance for detection of the SiO line features 


^max 

Teq = 2300 K 

Teq = 2500 K 

Teq = 3000 K 

4/itm 

- 

- 

340pc 

10/um 

30pc 

150pc 

430pc 

100/zm 

3.5pc 

8.7pc 

40pc 


5.3. Application to known HRSEs 

We show the simulated secondary-eclipse-depth spec¬ 
tra of four known close-in super-Earths, CoRoT-7b, 
Kepler-lOb, Kepler-78b, and 55 Cnc e in Fig. [T0l Within 
the lcr errors of their observed masses and radii, it is 
possible that those super-E arths are naked r ocky plan- 
ets. Note that, according to iDragomir et all • 12(11 1' i. the 
lcr upper limit on the measured 55 Cnc e’s density just 
reaches the pure rock regime. If so, they have mineral 
atmospheres, because of high T eq (see Table [3]). In the 
simulations, we have assumed 91 = 10 in 1-100 //rri and 
100 in 0.3-1 /jm. The values of the planetary and stellar 
parameters that we have used are listed in Table [3] 

The secondary eclipse depth becomes large, as R p /R * 
increases and T eq increases. Thus, the secondary eclipse 
depth of 55 Cnc e ( R p /R * = 2.12 Rq/Rq) is largest, 
while that of Kepler-lOb (R p /R* = 1.39 Rq/Rq) is 
smallest. Despite of its larger R p /R*, the secondary 
eclipse depth of CoRoT-7b is smaller than that of Keplcr- 
78b, because T eq of CoRoT-7b is smaller than that of 
Kepler-78b. 

Prominent absorption features are found around A = 
0.6 /im (Na), 0.77 /im (K), 4 /tm, 10 /tm and 100 /im 
(SiO). In particular, in the cases of Kepler-lOb and 
55 Cnc e, those features would be detectable by the near¬ 
future space-based observation that we have assumed in 
subsection 15.21 because the distances of these systems 
are shorter than (i?*/i? Q ) _1 (R p /2Rq) 2 d max from equa¬ 
tions o and C51) . 

iDemorv et all (|20l2l l have already reported the detec¬ 
tion of the secondary eclipse depth of 55 Cnc e with 


Table 3 

Assumed planetary and stellar parameters 



Corot-7b a 

Kepler-lOb 6 

Kepler- 78b c 

55 Cnc e d 

M p 

7.42 M© 

3.33 A/® 

1.69 M© 

8.09 M® 

R P 

1.58 R © 

1.47 Rq 

1.20 Rq 

1.99 i?0 

D 

0.0172AU 

0.0168 AU 

0.0092 AU 

0.0155 AU 

T eq 

2550 K 

3090 K 

3000 K 

2700 K 

R* 

0.82 Rq 

1.06 Rq 

0.74 Rq 

0.94 Rq 

T* 

5250 K 

5700 K 

5200 K 

5200 K 

distance 

150pc 

173pc 

- 

12.3pc 


0Howard et al. (2013), EPragomir ct al. (2014); Nelson ct al. 

f2014l ); lEhrenreich et all 120121 ) 


Notes. | a IHatzes et al. (2011); 

Bruntt et al. (2010); 

Leger et al .1 

(2009), ElDumusaue et al. 

(2014): Fressin et al 

(1201 U). 



Figure 10. The secondary eclipse depths theoretically predicted 
for (a) Corot-7 b, (b) Kepler-10 b, (c) Kepler-78 b and (d) 55 Cnc 
e, shown as a function of wavelength, assuming that they have the 
mineral atmosphere. The solid lines show the predicted secondary 
eclipse depth. The black dotted lines show the secondary eclipse 
depths of the blackbody, temperatures of which are indicated by 
”BB(T)”. We have assumed 91 = 100 in 0.3-1 /tm and 10 in 1-100 
/im and the planets’ parameters shown in Table [3] 

Spitzer Space Telescope at 4.5 /im. The measured bright¬ 
ness temperature of 2360 ± 300 K is consistent with our 
result shown in Fig. I10f d). although the error is large. 
Also, the SiO IR-absorption features of an HRSE that is 
brighter and closer than 55 Cnc e would be detectable 
with Spitzer in the IRAC 4.5 and 8 /tm channel. 

Using equations (l20l) and (l2ll> . we define the line con¬ 
trast to the total flux, aa) as 

Ch,(A,AA) = V\,A\ £A,AA- (24) 

The photometric accuracy should be smaller than 
C L , (a,aa) in order to detect the line spectra feature of the 
secondary eclipse depth. There are several proposed mis¬ 
sions that aim to characterize transiting exoplanets, such 
as Fast Infrared Exoplanet Spectroscopy Survey Explorer 
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Table 4 

Values of the line contrast 


the line contrast Cl(x,ax) 

Teq = 2300 K 

T eq = 2500 K 

Teq = 3000 K 

55 Cnc e 

Accuracy of EChO 

Na at 0.59 fim 

4.4ppm 

9.3ppm 

26ppm 

19ppm (9i = 100) 


- 

K at 0.77 fin 1 

6.4ppm 

14ppm 

36ppm 

26ppm (91 = 100) 


- 

SiO at 4 fim 

- 

- 

30ppm 

16ppm (9? = 50) 

12ppm (9? 

= 50) 

SiO at 10 fim 

4.6ppm 

24ppm 

67ppm 

54ppm (91 = 30) 

24ppm (9? 

= 30) 

SiO at 100 fim 

4.0ppm 

13ppm 

50ppm 

36ppm (9i = 30) 


- 


(FINESSE; e.g. Deroo jt all I 20 T 2 !) and t he Exop l anet 
Char acterization Observatory (EChO; e.g. iTinetti et al.l 
I2012T ). Here we compare our results with the accuracy of 
EChO. Figure fill shows the theoretical (solid line) and 
mock-EChO (bars) secondary-eclipse-depth spectra for 
a super-Earth with a mineral atmosphere, using 55 Cnc 
e’s parameters. Then, we create the mock-EChO spec¬ 
tra with 91 of 50 in 1-5 /un and 30 in 5-16 /im, and the 
photometric accuracy of ECh O in Chemical census mode 
(jWaldmann fe Pascalell2014l) Q. Our result demonstrates 
that EChO has enough accuracy to detect the SiO fea¬ 
ture of the mineral atmosphere of 55 Cnc e around 4 /im 
and 10 pm. 



Figure 11. The secondary eclipse depths theoretically predicted 
and the mock-EChO observational spectrum for 55 Cnc e with the 
mineral atmosphere. The secondary eclipse depth is shown as a 
function of wavelength. The solid line shows the predicted sec¬ 
ondary eclipse depth. The magenta bars show the mock-EChO 
observational spectrum dWaldmann &; Pascaldr2014D . We have as¬ 
sumed 91 = 50 in 1-5 fim and 30 in 5-16 fi m and 55 Cnc e’s pa¬ 
rameters: M p = 8.09 M 0 , R p = 1.99 R 0 , a = 0.0155AU, T eq 
2700 K #*=0.94 Rq and T*=5200K. 

Table 0] shows the line contrasts Cl of Na (0.59 fim), 
K (0.77 pm) and SiO (4 /im, 10 pm and 100 /im) for 
T eq = 2300 K, 2500 K, and 3000 K (91 = 100 for A < 
1 /iin and 91 = 10 for A > 1 /im), and 55 Cnc e’s pa¬ 
rameters (91 = 100 for A < 1 /im, 91 = 50 for 1 /im < 
A < 5 /im and 91 = 30 for A > 5 /im). The photometric 
ac curacy of EChO in Chemica l Census mode of 55 Cnc 
e (jWaldmann fc Pascaldf2014h is also shown in Table [4] 

3 IWaldmann K: Pascald (120141 ') showed the simulated EChO ob- 
servation of 55 Cnc e assuming a CO 2 and H 2 O atmosphere model. 
We use the noise level assumed in their Figure 9 because the depth 
of their spectra is similar to ours. 


Note that the values of Cl at 4 /im are not shown for 
T eq = 2300 K and 2500 K, because the Na, K and SiO 
features in those cases are indistinguishable with 91 = 10. 
We obtain that Cl,(\, aa) > 10 ppm for Na (0.59 /im), 
K (0.77 /im) and SiO (4 /im, 10 /im and 100 /im); for 
example, Cl,(\,a a) ~ 14 ppm at 0.77 /im and 24 ppm 
at 10 /im for T eq = 2500 K and Cl,(x,a a) ~ 26 ppm at 
0.59 /im, 34 ppm at 0.77 /mi, and 67 ppm at 10 /im for 
T eq = 3000 K. The required photometric accuracy to de¬ 
tect the SiO features at 4 and 10 /im is achievable with 
the Chemical Census mode of EChO for 55 Cnc e. 

5.4. How many HRSEs are detectable via near-future 
observation? 

The Kepler space telescope discovered 4250 planet can¬ 
didates (http://exoplanets.org, as of 30 Nov. 2014). The 
Kepler Objects of Interest (KOIs) include super-Earth- 
size objects with short periods (i.e., HRSEs), in particu¬ 
lar, 751 objects with R p < 2 Rq orbiting stars with T e © 
= 5300-6000 K. Among them, there are 60 KOIs with 
T eq > 2300 K. If they are naked rocky planets, their min¬ 
eral atmospheres can be identified by secondary-eclipse 
observations with sufficient signal-to-noise ratios. We es¬ 
timate the expected number of such HRSEs detectable 
by near future observations. 

Assuming a uniform distributio n of G- type stars, the 
number of the detectable HRSEs, Nurse, is given by 

_ 47 r 2 

Ahrse ~ rNc— [dmax {(iS/IV)aa}] j (25) 

where r is the occurrence rate of HRSEs with T eq > 
2300 K hosted by G-type stars and Nq is the number 
density of G-type stars. The HRSEs with T eq > 2300 K, 
which should exhibit the prominent emission features 
(Fig. 0), correspond to pl anets with peri od P < 2 days 
around G-type stars. iFressin et al.1 (1201311 estimated, us¬ 
ing KOI host stars, the occurrence rate of planets with 
1.251?© < R p < 2i?© hosted by G-type stars to be close 
to 0.35. They also estimated the occurrence rate of plan¬ 
ets with 1.251?© < Rp < 2i?© and 0.8 < P < 2 days 
to be 1.7(±0.3) x 10 . Using their results, the occur¬ 

rence rate of the planets with 1.251?© < R p < 21?© and 
0.8 < P < 2 days hosted by G-type stars is estimated to 
be r ~ 6 x 10 -4 . Here we use r = 6 x 10~ 4 as a fiducial 
value of the occurrence rate. 

In section El we estimated 

dmax{(S'/A) A A = 3 /t/a,aa} to be 430 pc (T eq = 
3000 K), 150 pc (T eq = 2500 K) and 30 pc (T eq = 
2300 K) for 10 /mi observation (see Fig. [9] and Ta¬ 
ble 0. These HRSEs with T eq = 2300 K, 2500 K 
and 3000 K have the orbital period P ~ 0.93, 1.6 
and 2.0 days, respectively. Hence, we take the sim¬ 
ple mean [(430 3 + 150 3 + 30 3 )/3] pc ~ 300pc as 
a fiducial value of the limiting distance for HRSEs 
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with 0.8 < P < 2 days. Inputting r = 6 x 10 -4 
% and Ng ~ 6xl0~ 3 pc~ 3 (see a recent review by 
iTraub fc Onncnheimerl l2Qllh in equation (|25]>. we find 
that within 300 pc there are approximately 400 HRSEs 
whose 10 fim SiO features are detectable . Taking the 
transiting probability of ~ 0.2 (~ Rq/D for P = 1.3 
days) into account, one will be able to detect the 10 /im 
SiO features from ~ 80 transiting HRSEs by a 5 m 
space telescope equipped with a photon-noise limited 
instrument. 

6. MODEL VALIDATION AND FUTURE STUDY 

6 .1. Scattering 

In this study, we have ignored the effect of scatter¬ 
ing on the thermal structure of the mineral atmosphere. 
There are two possible sources for scattering. One is 
Rayleigh scattering by the major gases. The cross section 
for Rayleigh-scattering by gas species A, ignoring 

the p olarizability anisotropy , are approximately given as 
(e.g., ICRC Handbookll201l(i 

= ^(2t Tvfa\, (26) 

where a a is the static average el ectric dipole polarizabi l- 
ity of gas species A. According to lCRC Handbook! (|201 11 1. 
a a for ground-state atoms and diatomic molecules in the 
unit of 10 -24 cm 3 are 24.11 for Na, 43.4 for K, 0.802 for 
O, and 1.5689 for O 2 . While oa for SiO is unavailable in 
the literature, to our knowledge, it should not differ by 
several orders of magnitude from those of the other gases. 
Comparing the cross sections for Rayleigh scattering by 
those gases with the total absorption cross section calcu¬ 
lated in section[3l we have confirmed that the formers are 
smaller by several orders of magnitude than the latter in 
most of the wavelength region of interest, except for the 
narrow ranges between 0.1 and 0.15 /im and between 0.3 
and 0.45 ^m. We have integrated the thermal structure 
incorporating the effect of the Rayleigh scattering by Na, 
which results in only a few percent difference in temper¬ 
ature. Also, the Rayleigh and Mie scattering by cloud 
particles can affect the thermal structure of the mineral 
atmosphere, if present. However, since the atmospheric 
temperature is always higher than the vapor pressure of 
rock, as shown in Fig.|4l no condensation cloud is formed 
in the atmosphere. 

6.2. Non-LTE effect 

The approximation of local thermal equilibrium (LTE) 
is invalid at low pressures. In this study, we have con¬ 
sidered a pressure range of > 1 x 10 -8 bar. For LTE 
to be attained, thermalization due to atomic/molecular 
collisions must occur faster than radiative cooling. The 
collisional timescale, t co ii, is ~ kT /(Pi><t co ii), where a co u 
and v are the collisional cross section and thermal veloc¬ 
ity of molecules. For SiO, since <t co ii ~ 3 x 10 -15 cm 2 , T 
~ 3000 K, v « 8 x 10 4 cm/s and the typical radiative- 
cooling timescale (i.e., the typical lifeti me o f excited 
states) t coo i « 1 x 10 -2 s dBarton et al.l I2013h . LTE is 
valid for P > 2 x 10 -7 bar. Since the emission from 
the low pressure region of P < 2 x 10 - ' bar makes little 
contribution to the secondary eclipse depth, the LTE ap¬ 
proximation is valid for the SiO-dominated atmosphere 
of special interest in this study. 


6.3. Photodissociation of SiO 


Our equilibrium calculations show that the major gas 
species in the mineral atmosphere are Na, K, Fe, O, O 2 , 
Si and SiO (section U. Although having been ignored 
in this study, photodissociation of SiO would occur in 
such an atmosp here by irradiation of stellar UV (e.g., 
lVorvpacvlll98lD. 

iJolicard et al l (119971) provided the photodissociation 
cross section of SiO, <7di s , in the wavelength range be¬ 
tween 0.153 /im and 0.157 /im, where the cross section 
decreases monotonically with increasing wavelength from 
~10 -17 cm 2 to 10 -19 cm 2 . The energy of UV rays of 
A > 0.157 /im is not high enough to dissociate SiO sig¬ 
nificantly; instead, the electron-transition of SiO is dom¬ 
inant, as discussed in section [3] Also, for A = 0.152- 
0.163 /im, the electron-transition cross section of Si is 
larger than 10” 17 cm 2 and thus larger than <Tdi s . There¬ 
fore, it is the UV ray of A < 0.152 /jm that dissociates 
SiO efficiently in the mineral atmosphere. Although no 
other published data for <Tdi s is available, <Tdi s is unlikely 
to decrease with decreasing wavelength. Below, we make 
a rough estimate on the impact of SiO photodissociation 
on the atmospheric structure and the emission features. 

Using Udis > 1 x 10 -17 cm 2 with xsio = lx 10 -3 , 
Ji = 20 and g = 2500 cm/s -2 , we calculate the SiO- 
photodissociation optical depth, which is defined by 


^dis 


^SiO^dis p 

Jimug 


(27) 


where toh is the proton mass. Then, we find that Tdi s = 1 
at P < 8 x 10 -6 bar. This means that the atmospheric 
structure for P<8xl0 -6 bar is affected significantly 
by photodissociation of SiO. Since the emission features 
shown in Fig.[]]arise from the region of P > 3 x 10 -5 bar 
(see Fig. E}, the emission spectra are unlikely to be af¬ 
fected by the photodissociation. More detailed exami¬ 
nation requires photo-chemical calculations done consis¬ 
tently with radiative transfer calculations, which must 
be a future study. 


6.4. Tidal heating 

Our atmospheric models do not include the effect of 
tidal heating, discribed in se ction [4] (Fn = 10 W/ m 2 be¬ 
ing assumed). According to iHenning et al.l ( 2009 ). tidal 
heating dominates stellar irradiation under somewhat ex¬ 
treme circumstances such that the planetary eccentricity, 
e, is larger th an 0.1 and the orbi tal period is less than 
2 days. Also, iMiguel et al.l (120111) showed that addition 
of tidal heating hardly changes the surface temperature 
of HRSEs for e = 0.01 and 0.001 and the quality fac¬ 
tor of the planet, Q , of 200 and 20000. Indeed, although 
not shown, we have confirmed that there is no significant 
change in the atmospheric structure for T eq = 3000 K, 
by calculating the temperature profiles with F 0 = 10 3 - 
10' W/m 2 , being equivalent to tidal heat flux for e = 
0.01 and 0.001 and Q = 200 and 20000. 

6.5. Temperature jump at the bottom 

Allowing the temperature jump would affect the com¬ 
position of the atmosphere. If we included the thermal 
conduction, we would have no temperature jump at the 
bottom in the atmosphere, and would have higher T g 
and a thicker atmosphere than our model in this study, 
which is more detectable. The HRSE atmosphere is thin 
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enough for radiation to heat the atmosphere more ef¬ 
fectively than thermal conduction. Thus, although the 
thermal conduction would change T g and thus the com¬ 
position of the atmosphere, it would be ineffective in 
changing the emission spectra of HRSEs. 

6 .6. Impact of potassium on thermal structure and 
detectability 

As shown in section [2j our chemical-equilibrium cal- 
culations yield much m ore abundant potassium than 
iSchaefer fc Feelevl (12009H did. Although not shown, we 
have confirmed that such difference in K abundance has 
only a small impact on the thermal structure of the 
atmosphere: Our calculations have demonstrated that 
the maximum differences in temperature are ~2.86%, 
~2.27% and ~1.89% for T eq = 3000 K, 2500 K and 
2300 K, respectively. Thus, the secondary eclipse depths 
in the IR region are hardly affected. 

A direct impact on the secondary eclipse depth is found 
at A = 0.77 /im. Detecting potassium may be useful for 
constraining the magma composition, because the potas¬ 
sium abundance is different by one order of magnitude 
from MORB to the upper crust, as shown in Figure [2] 

6.7. Other observation methods to detect mineral 
atmospheres 

Other observation methods for atmospheres of close-in 
exoplanets include transmission spectroscopy and plane¬ 
tary radial velocity (PRV) measurement, in addition to 
secondary eclipse observation. It would be, however, dif¬ 
ficult to detect the HRSEs’ atmospheres via transmis¬ 
sion spectroscopy of the planetary disk because of their 
small scale height. The part of the HRSE’s atmosphere 
observed by transmission spectroscopy is the rim where 
horizontal heat and composition transport between the 
day side and night side would be important. The ex¬ 
amination requires 3-D calculations of atmospheric cir¬ 
culation done consistently with radiative transfer calcu¬ 
lations, which must be a future study. 

Recently, several molecules in the dayside spectra 


of hot Jupiters have been successfully 

detected us- 

ing the high dispersion spectrograph (Brogi et al. 20121 

Rodler et al.1 12013: Ide Kok et al.l |2013t 

Rodler et al.l 

20131 Birkbv et al.1120131 iLockwood et al.1 

20llT This 


method utilizes the cross-correlation analysis of the spec¬ 
tra and the model template and directly detects PRV. 
Here, we estimate the feasibility of the PRV detection of 
the mineral atmosphere of HRSEs. For the most opti¬ 
mistic case, the signal-to-noise ratio of the PRV is deter¬ 
mined by the photon noise, 

{S/N) p RV = \ZNm^C(S/N) A xC L ,(x,A\), (28) 

where Mine is the effective number of molecular lines 
available for the detection and ( is th e total efficiency 
of th e high-dispersion instrument (see iKawaha ra et al.1 
I2014T) . Although the determination of the precise value 
of Mine requires full simulations of the spectra includ¬ 
ing the terrestrial transmission and the night airglow, we 
here perform 1st order-of-magnitude estimate of the de¬ 
tectability of the Na and K lines here, assuming a future 
30 m-class ground-based telescope. 

The number of Na and K lines within the secondary 
eclipse depth spectrum of an HRSE for T eq > 2300K, 


Vii ne , is ~ 100, the line contrast, C l jx,a\), is ~ 0-5 
and the mean value of those depths, ex is ~ 10 -4 , in 
the IR wavelength region (1-10/iin), as shown in Fig. [6] 
The resolving power of the instrument, 91, is required 
to be large enough to resolve those Na and K lines. To 
know this value of 91, we have reintegrated the secondary 
eclipse depth for T eq = 2500K around only several Na and 
K line peaks in the IR wavelength region (3.70-3.75/im) 
by line-by-line calculation. We have confirmed that 91 = 
10 5 is enough to resolve the Na and K lines. Thus, we use 
91 = 10 5 below as a fiducial value of 91. From equation 
K5l) . the distance of the system within which at least 
one transiting HRSE exists is ~ 70pc. Assuming a 30 m 
telescope with ( = 0.1, we obtain the expected S/N for 
detecting the Na and K lines, 


(S/N) P ny ~ 
, Mine C 


ex 70pc i?* Cl,(x, aa) 
~d~lR ,q R5 

t ohs 10 5 7t\B(T*)//i 


io - 4 

A 


100 0.1 7rl5 2 m 2 lOhr 91 10 26 m _2 s _1 


(29) 


where 10 26 m~ 2 s _1 corresponds to the mean value of 
7 rB(T*)/h in the wavelength region between 1 and 10/rni, 
assuming T* = 6000 K. This order-of-magnitude esti¬ 
mate suggests that it would be possible to detect Na and 
K in the HRSEs’ atmospheres via PRV measurement in 
the very large telescope era. 


7. SUMMARY AND CONCLUSIONS 

In this study, we have examined the detectability of 
hot rocky super-Earths that have atmospheres in the 
gas/melt equilibrium with the underlying magma oceans 
via secondary-eclipse observations. First, we have done 
gas-melt equilibrium calculations and showed that the 
major atmospheric species are Na, K, Fe, SiO, O, and 
O 2 (§0). We have called this atmosphere the mineral 
atmosphere in this study. 

Next, for those species, we have compiled the absorp¬ 
tion line data available in literature, and calculated their 
absorption opacities in the wavelength region of 0.1- 
100 pm (§[31). We have found that the SiO absorption 
dominates in the UV and IR wavelength regions, while 
strong absorption by Na and K is found in the visible 
region. In this atmospheric gas, the opacity is larger in 
the UV and visible regions than in the IR region. 

Third, using those opacities, we have integrated the 
thermal structure of the atmosphere (§ IH). The most 
important finding is that thermal inversion occurs in the 
atmosphere with relatively high equilibrium temperature 
(i.e., pj T eq > 1900 K), because the UV and visible ab¬ 
sorption by SiO dominates the IR absorption by SiO and 
others. The presence of thermal inversion has a crucial 
impact on the detectability via secondary eclipse obser¬ 
vation. 

Finally, by calculating the ratio of the planetary and 
stellar emission during secondary eclipse, we have quan¬ 
tified the detectability (§ ©. We have found prominent 
emission features due to SiO at 4, 10, and 100 pm. Com¬ 
paring our predicted spectrum with a mock spectrum 
we made, we have demonstrated that the mineral atmo¬ 
sphere of a hot rocky super-Earth of radius 2 with 
T eq > 2300 K orbiting a solar analog would be detectable 
by a 5-m space telescope. 
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The Transiting Exoplanet Survey Satellite (TESS) is 
a planned space telescope for NASA’s Explorer program 
and is dedicated to an all-sky transit survey of exoplan¬ 
ets (iRickcr ct alJl2010h . It is predicted that TESS will 
discover more than 100 transiting Earth-sized exoplan¬ 
ets with short orbital periods. After discovery of nearby 
HRSEs, follow-up observations by satellites with large 
telescopes such as JWST, Spitzer, FINESSE and EChO 
will open new exotic worlds of the mineral atmospheres. 
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APPENDIX 

APPENDIX A: CALCULATION METHOD OF THE MEAN OPACITIES 

In this study, we use mean opacities for simulating the atmospheric structure and emission spectra, instead of line- 
by-line calculations. While we integrate the line opacities of O, O 2 , Na, K and Si with wavelength resolution high 
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enough to reproduce line profiles, we use low resolution grids to calculate the mean opacity of SiO, because there are 
too many absorption lines to reproduce the line profile of SiO in a practical CPU time. 

Figure [T2] presents how to calculate the mean opacity of SiO when an interval between adjacent lines is less than 
1 /im. When an interval between adjacent lines is smaller than 2Aud and 2(r r +Tw) (see the yellow-colored region 
in Figure IT2l) . we divide the interval equally into quarters. On the other hand, when an interval is larger than 2Avd 
or 2(r r + Tw) (see the blue-colored region in Figure [12]), we select a total of five grids: vo,i, r'o.i+i, the larger of 
vo,i + 2,Avd and iso ,i + 2(r r + rw) (= vi), the smaller of vqj+i — 2Avd and vo,i+i — 2(r r q-Tw) (= v r )i and a midpoint 
between and t'o.i+i (— v m)i where uq^ is the central frequency of fth spectral line. When we integrate the line 
opacities around the line center, we use a trapezoidal integration method. On the other hand, when integrating line 
opacities between U( (i/ r ) and u m , we approximate the integration as the sum of the areas of the shaded boxes. 

When an interval is larger than 1 /im, we use the line opacities at u o i , v m and 50 grids between those points; those 
grids are selected so that the size of each interval logarithmically increases with frequency. The integration is done in 
the same way above. 

When simulating the atmospheric thermal structure, we adopt the harmonic mean opacities. This is because the left 
side of equation (flOl) contains the Planck’s function B V (T ) and its right side contains the opacity in the denominator of 
dF±/dr u . On the other hand, when calculating the planetary emission spectra, we use the arithmetic mean opacities. 
This is because equation (1181) includes the Planck’s function multiplied by the exponential of the negative opacity and 
the opacity (i.e., B\e~ Tx dT\), which can be expressed as the sum of B\k ™ (m = 1, 2, 3, • • •). 



Figure 12. Schematic line profile of SiO (red line) with illustration for selecting grids (red points) to calculate the SiO mean opacity. In 
the yellow-colored (blue-colored) region, the interval between two adjacent lines is smaller (larger) than 2A vp and (or) 2(IY + T^y). For 
example, the opacity integrated between and ^o,i+i is approximated to be the sum of the areas of three shaded boxes. 















